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1 Introduction 


Micromechanics techniques can be employed to model the individual constituents 
within the composite. Typically, a repeating unit cell (RUC) in the composite mi- 
crostructure is identified and analysis is performed on that RUC assuming periodic 
boundary conditions. The response of a point in a continuum is determined assum- 
ing an infinite array of the RUCs. However, representative volume element (RVE) 
methodologies exist which incorporate applying non-periodic boundary conditions to 
a subvolume that accurately represents the composite microstructure. The RVE is 
meant to represent the actual microstructure of the continuum, and the size of the 
features of the microstructure is preserved. Micromechanics can be utilized to provide 
the homogenized composite stiffness, or they can be used to model damage and fail- 
ure within the constituents. If utilized for the latter, the global mechanisms can arise 
through the natural evolution and interaction of the mechanisms in the constituents 
in the micromechanics model. Numerous micromechanical frameworks exist that en- 
compass analytical, semi- analytical, and numerical techniques. An expansive review 
of many micromechanics theories is given in Ref. [Aboudi et al. (2012)]. In this work 
a discretization insensitive continuum damage model (CDM) is implemented within 
the high-fidelity generalized method of cells (HFGMC) to model transverse cracking 
and compressive shear banding in a unidirectional fiber-reinforced polymer matrix 
composite (PMC). 

The method of cells (MOC) developed by Aboudi (1991) discretized a rectangu- 
lar composite RUC into four subvolumes, called subcells. One of the subcells was 
occupied by the fiber material and the rest were occupied by the matrix. Linear 
displacement fields were assumed in each of the subcells. Displacement and traction 
continuity conditions were enforced, in an average integral sense, at the subcell in- 
terfaces, along with periodic boundary conditions at the RUC boundaries to derive 
a set of equations that would yield a strain concentration matrix that could, in turn, 
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be used to obtain the local subcell strains from the applied fields. Following deter- 
mination of the subcell strains, the subcell stress are readily calculated using the 
local constitutive laws, and volume averaging can be used to obtain the homogenized 
thermomechanical properties of the composite. MOC was later extended to the gen- 
eralized method of cells (GMC) by Paley and Aboudi (1992) which accommodated 
any number of subcells and constituents in two-periodic directions. Aboudi (1995) 
adapted the formulation to accommodate triply-periodic materials. Finally, Aboudi 
et al. (2001) developed HFGMC which utilized second order displacement field ap- 
proximations in the subcells, rather than linear. Haj-Ali and Aboudi (2009) showed 
the local elastic field accuracy produced by HFGMC corresponded very well to FEM; 
whereas, Aboudi et al. (2003) compared the accuracy of the fields in inelastic phases 
to composite cylinder assemblage (CCA). Bednarcyk et al. (2004) utilized HFGMC to 
model fiber-matrix debonding in metal matrix composites (MMCs), and Bednarcyk 
et al. (2010) implemented a multi-axial damage model in HFGMC. Reformulations, of 
GMC and HFGMC, which reduced the total number of unknowns in the problem were 
introduced by Pindera and Bednarcyk (1997) and Aboudi et al. (2012), respectively. 

The generality of the GMC and HFGMC formulations admit any constitutive 
behavior at the subcell level. However if the response of the subcell material exhibits 
post-peak strain softening the tangent stiffness tensor of the subcell loses positive 
definiteness. This leads to pathologically mesh dependent behavior. 

Bazant and Oh (1983) developed the smeared crack band approach which intro- 
duced a characteristic element length into the post-peak softening damage evolution 
formulation. The tangent slope of the softening stress-strain curve was scaled by 
the characteristic length to ensure that total strain energy release rate (SERR) upon 
complete failure (i.e. zero stress) is always equal to the prescribed fracture tough- 
ness, regardless of the element size. In the original formulation, the band was always 
oriented perpendicular to the direction of maximum principal stress; thus, the crack 
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band always advanced under pure mode I. de Borst and Nauta (1985); Rots and 
de Borst (1987) later reformulated the model to incorporate a fixed crack band that 
evolved under mixed-mode conditions. Both formulations employ linear degradation 
schemes. Later, Camanho et al. (2007) incorporated more sophisticated initiation 
criteria to predict the onset of mixed-mode crack bands. All of these smeared crack 
formulations assume linear elastic behavior up to the initiation of the crack band, 
followed by immediate post-peak strain softening. However, Spencer (2002) coupled 
pre-peak plasticity with crack band post-peak strain softening in modeling failure of 
concrete. 

A variation of the crack band model developed by Bazant and Oh (1983) for con- 
crete structures is implemented here within the HFGMC micromechanics framework 
to model mesh objective failure of continuous fiber-reinforced polymer matrix com- 
posites. The focus is restricted to the microscale to evaluate the capabilities of the 
smeared crack band model to predict progressive failure evolution within a composite 
microstructure. Thus, a detailed, multiple-fiber RUC is chosen demonstrate these 
capabilities. HFGMC is utilized as the micromechanics platform because the level of 
fidelity required to accurately model failure evolution in a complex microstructure is 
not provided by GMC. 

In Section 2 the formulation for HFGMC is given, followed by the implementation 
of the smeared crack band model in Section 3. Two scenarios are considered to deter- 
mine the mode in which the cracks within the crack band grow. If the principal stress 
that has the highest magnitude is tensile, it is assumed that it is more energetically 
favorable for the crack band to form perpendicular to the maximum principal stress 
and for the cracks within the band to advance under mode I conditions. Conversely, if 
the the magnitude of a compressive principal stress is higher than the other principal 
stresses, the cracks within the crack band evolve under mode II conditions (due to 
internal, Mohr-Coulomb friction) and are oriented with the plane of maximum shear 
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stress. In Section 4 the objectivity of the post-peak softening behavior of the model 
with respect to the subcell mesh is shown. An example is presented in Section 5 
that consists of an RUC composed of 13, randomly placed fibers that is subjected to 
transverse tension and compression. The model is verified against experimental data 
in Section 6. 

The motivation behind the following developments is to arrive at a failure method- 
ology that is capable of accurately capturing localization in a multiscale model without 
any spurious dependencies. This paper is intended to present results that verify a 
model for such an application. In Section 7 a multiscale methodology is introduced 
that ensures mesh objectivity and addresses the issues presented by Bazant (2007). 

2 The High-Fidelity Generalized Method of Cells 

HFGMC, first introduce by Aboudi et al. (1999) resolves some the shortcomings of 
the original GMC. The biggest of which is the lack of normal-shear coupling in GMC 
which results from the enforcement of displacement and traction continuity in an 
average, integral sense and a linear approximation of the displacement fields in each 
subcell. HFGMC employs quadratic displacement approximations. However, this 
requires more equations than can be produced using displacement and traction con- 
tinuity. Thus, the zeroth, first, and second moments of equilibrium are used to solve 
the problem. 

It is assumed that a composite microstructure can be represented as a collection of 
triply periodic RUCs containing a general number of constituents, as shown in Figure 
1. The RUC is then discretized into N a x Np x N 1 rectangular, prismatic subcells, as 
exhibited in Figure 2. Each of these subcells are occupied by one of the constituents in 
the composite. The number of subcells and the materials occupying each subcell are 
completely general. For a two-phase fibrous composite any desired micro-architecture 
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can be represented by occupying each subcell with either a matrix or fiber constituent. 
Since the microstructure does not vary along the axial, fiber-direction, a unidirectional 
composite can be modeled using a doubly-periodic RUC (Figure 3), where the X 2 -X 3 
plane is discretized. Herein, doubly-periodic formulations of HFGMC [ Aboudi et al. 
(1999)] will be employed; the reader is referred to Ref. [ Aboudi et al. (2012)] for the 
fully 3D formulations, as well as a more efficient reformulation. 

2.1 Governing Equations 

The 3D equilibrium equations for subcell by (see Figure 4) in a doubly-periodic RUC 
(Figure 3) are given by 


sCT + ftdf 1 + a 3 4 w = 0 


(i) 


Since there is no variation in the x , -direction for a doubly-periodic material, 0\ = 
d/dxi = 0. Therefore Equation (1) becomes 


&d? 7) + s,d? 7) = o 


Averaging Equation (2) over the volume of the subcell yields 


jC7) _|_ tA/A) _ Q 
J 2i(00) ^ -“-31(00) — U 
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The surface tractions are given by 


(2u±(/3 7 ) = _L [ v 2 ( '$» = ±h\ d ^) 


i i 2i 

^7 J —l-y/2 


(3) t ±m = J_ 


rh " /2 ( X { : ] = )dxi? ) 


h P J-hp / 2 3 * 


(6) 


(7) 


where “(2)” and “(3)” are superscripts that indicate the applicable normal direction 
to the subcell surface and “±” indicates the sign of the normal. Similarly, the first 
moment of equilibrium can also be averaged over the subcell volume resulting in 


j(Pl) _ q(Pl) _ n 

J 2i(l0) °2i(00) U 


(8) 


and 


zy(^7) _ q(Pl) _ r, 

-^31(01) °3i(00) — U 


where , is an average stress quantity, and 
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(10) 

( 11 ) 

( 12 ) 


ij{mn) hpL, J_ h0/2 J_ h/2 

Finally, the second moment of equilibrium can be averaged over the subcell volume, 
then utilizing integration by parts leads to 


hi 


\ 1 j(Pl) I ^ P Tf(Pl) 
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_ n _ O C (P'1) - 0 
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(13) 
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Substituting Equations (13) and (14) into Equation (3) gives 


j(07) q(Pl) 

J 2i(00) 7,2 °2i(10) 

,l P 


id: 


(^7) _ 


12 


Q\r 

LJ q • 


(/? t ) 


3i(00) p ^3i(01) 


Utilizing Equations (15) and (16) in Equation (14): 


1 0(^7) 1 J_ 0(^7) 

7,2 °2i(10) _r 72 J 3i(01) 
,< 73 s 


0 


(15) 

(16) 


(17) 


Equation (17) represents the average form, in HFGMC, of the equilibrium equations, 
Equation (2), of subcell /3y within a doubly-periodic composite RUC. 


2.2 Second Order Displacement Approximation 

Contrary to GMC, HFGMC employs a quadratic expansion of the subcell displace- 
ments 


u, 


(07) 


— t-jjXj 


| Wjgg | + x^’W) 


^ ( 3x{ 0) r- 1 W, 


h 2 a 


7 ( 10 ) 

r(07) 
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7 ( 01 ) 


1 
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7 (^ 7 ) 
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where e, 3 are, as in GMC, the average global strain components, and are 

microvariables describing the spatial variation of the displacement field and are de- 
termined by enforcing interfacial and periodic conditions. 


2.3 Constitutive Relations 


The local strains in subcell /I7 are related to the global strains and local displacement 
fields through 


(07) 

ij 


tj 


+ 1 (Vf 7) + d.uT) 
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Employing the local constitutive law, 


AM) _ MM) AM) 

a ij ~ ° ijkl t ki 


( 20 ) 


ignoring plastic and thermal effects (for HFGMC formulation with plasticity and 
thermal strains please see Ref. [ Aboudi et al. (2003)]) in Equation (12) along with 
Equations (18) and (19), the relationship between zeroth order stress components 
and the microvariables is determined to be: 

4i(oo) = + C[P (w^ + e 22 ) + (w^fo 7 !) + e 33 ) 

+ ^ 7) (<o?) + <4 + 2e 23 ) + C ™ > (W^ + 2e 13 ) (21) 
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^55(00) = + c ' 25 7) (^ifio) + ^22) + Cg 7) (w^l) + £33) 

+ c[l l] (V^ + Wg^o) + 2 623^ + cg 7) (hh^Ol) + 2 ei 3 ^ 
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where the contracted notation for the stiffness tensor components C'P 7 ' 1 is used, 
ilarly the first stress moments are 
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2.4 Displacement Continuity 

The unknown microvariables (the volume- averaged displacement vector Wq/oo) and the 
higher order terms iu7i)) can be calculated from the governing equation, Equation 
(17), subcell interfacial continuity conditions, and periodicity conditions. 

The periodic boundary conditions for an RUC are given by 


Ui \x 2=0 — U i\x 2 =H 


(37) 


(J ‘Zi\x 2 = 0 — i\x 2 =H 


U i\x 3 = 0 — U i\x 3 =L 


| X3 _q 0's i | X3 


(38) 


(39) 


(40) 


The displacement continuity conditions are enforced, as in GMC, in an average (in- 
tegral) sense, yielding: 
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u. 


m 
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(41) 


dx[ f \ 0 = 1 ,...,Np 


(42) 


Using the definition of the subcell displacements, Equation (18), in the average dis- 
placement continuity conditions, Equations (41) and (42), gives: 


wdU _ _ MwdU _ w v*m) _ HwW _ % w(^7) 1 N 
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Equations (43) and (44) represent the displacement continuity conditions necessary 
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to enforce periodicity of the RUC. Similar relations can be derived to enforce dis- 
placement continuity across adjacent subcell boundaries, 


ti /(/U) I ^P \\r(P*t) | Tjb/U) TJ/O^+U) | ^/3+l Ti/ft’+l'/) | Ti/t^+U) 

^i(OO) ■+■ 2 ^*(10) + 4 ” i(20) — ^i(OO) + 2 ^*(10) + 4 ” i(20) ’ 
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(45) 
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^i(oo) + 2 *( 01 ) + 4 h02) “ hoo) + 2 *( 01 ) + 4 *(° 2 ) ’ ( 46 ) 

/5 = 1, . . . , IVg, 7 = 1,..., Nj — 1 

Note that Equations (43)-(44) can be relaxed to incorporate the effects of local, 
interfacial debonding [ Bednarcyk et al. (2004)]. 


H 


l. 


H 


2.5 Traction Continuity 

Similarly using Equations (6) and (7), the average implementation of the stress peri- 
odicity conditions (Equations (38) and (40)) is 
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Using Equations (4) and (5) with Equations (10) and (11) gives 
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Substituting Equations (8), (9), (15), and (16) 
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{3) tf m = slim ± r/mi)’ P = h...,N s , 7 = 1, . . . ,TV 7 (52) 

Thus, the average stress periodicity conditions, Equations (47) and (48), become 
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Equations (53) and (54) represent the stress continuity conditions at the boundaries 
of the RUC necessary to enforce periodicity. Similarly, traction continuity between 
adjacent subcells is given by: 


q(P 7) q(P 7) c(/ 3 + 1 7) 

°2*(00) h ° 2i(10) °2i(00) h 
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oi(00) 1 ^3i(01) 
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(55) 
. 

(56) 


2.6 Summary of Global Equations 

The number of unknown microvariables in the displacement expansion, Equation (18), 
is 15 NpN~,. The governing equations, Equation (17) provide 3 NpNy relations for the 
unknown microvariables while Equations (43) and (44) provide 3(Np + AT,) relations, 
Equations (45) and (46) provide 3[(-/Vg — l)iV 7 + Np(Nj — 1)] relations, Equations 
(53) and (54) provide 3 (Np + N 7 ) relations, and Equations (55) and (56) provide 
3 [(Np — 1 ) Nj + Np(Nry — 1)] relations, which yields a total of lSNpN^ linear equations 
for the unknown microvariables. 

This system of equations can be expressed as 

KU = f (57) 
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where K, the structural stiffness matrix, contains geometric and mechanical consti- 
tutive information from each of the subcells in the RUC. The displacement vector U 
contains the unknown displacement microvariables for each subcell. 


u = [u ( 11 ) ,...,u (N » N -' ) ; 


(58) 


and 

Tj07) - [w(W wOO w (/?7) W (/?7) W (/?7) l (Aql 

U — L VV (00)’ VV (10)’ VV (01) > VV (20)’ VV (02) J l 0y l 

where the vectors contain the three components of that displacement mi- 

crovariable for that particular subcell 


W 


(Pi) __ 

(mn) 


' W (07) W (07) W (Pl) ' 

VV l(mn)’ VV 2 (mn)’ 3 (mn) 


(60) 


The force vector f contains details of the applied average strains e tJ . To prevent 
rigid body motion, the displacement of the corners of the RUC are pinned and those 
equations are eliminated from Equation (57). 

Solving Equation (57) yields a strain concentration matrix that relates the local, 
average subcell strains to the global, applied, average strains. 


e «J-r> 



( 61 ) 


Once the local strains are obtained, the local stresses are readily determined through 
the local constitutive laws, Equation (20). The average global stresses are simply the 
volume average of the local stresses. 


cr 


1 

Til 


Np 

V7 ^ (/37) 


/3=1 7=1 


(62) 


where <r^ 7 ^ are the average subcell stresses. 
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The global, composite constitutive relationship is taken to be 


v = C*e (63) 

Thus, the effective elastic stiffness matrix for the composite is given by substituting 
Equations (61) and (62) into Equation (63). 

Np Aq 

7 C^A^ } (64) 

/ 3=1 7=1 

It should be noted that HFGMC was reformulated using the displacement continuity 
conditions to reduce the overall number of unknown microvariables, thus reducing 
computer memory requirements and computational cost. The reader is referred to 
[Bansal and Pindera (2004); Aboudi et al. (2012)] for details on this reformulation. 

3 Modeling Constituent-Level Post-Peak Strain 
Softening with the Smeared Crack Band Ap- 
proach 

HFGMC is an efficient (relative to fully numerical methods), useful tool for modeling 
details of the microstructure of a composite material. Additionally, it is readily 
amenable for implementation into a multiscale framework. Although, physics-based, 
discretization objective, progressive failure constitutive models must be in place for 
the constituents of the composite to accurately predict the response of a structure 
that is a damaging. For pre-peak loading (i.e. positive-definite tangent stiffness 
tensor), there are a multitude of non-linear elasticity, plasticity, continuum damage 
mechanics, and visoelastic/plastic theories available that can predict the evolution of 
the appropriate mechanisms in the composite. However when the local fields enter 
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the post-peak regime of the stress-strain laws, most of these theories breakdown in 
a numerical setting and display pathological mesh dependence [ Bazant and Cedolin 
(1979); Pietruszczak and Mroz (1981)]. 

The lack of positive definiteness of the elastic, or inelastic, tangent stiffness tensor 
leads to imaginary wave speeds in the material. The longitudinal wave speed in an 
isotropic material is given by 


c L 


E( 1 - v) 
p( 1 + /')( 1 - 2^) 


(65) 


where Cl is the wave velocity, E is the Young’s modulus of the material, v is the Pois- 
son’s ratio, and p is the material density. A one- dimensional approximation yields 
v = The existence of an imaginary wave speed results in a boundary value 

problem that is ill-posed Bazant and Cedolin (1979); Pietruszczak and Mroz (1981); 
deBorst (1987); Bazant and Cedolin (1991). Physically, a material must posses a 
positive-definite tangent stiffness tensor, and in fact, at the micro-scale the material 
tangent stiffness tensor always remains positive-definite. However for practical pur- 
poses, engineers must model structures at scales much larger than the flaws in the 
material, and the homogenized continuum representation of a material with nucleation 
and propagation of discontinuities, such as cracks or voids, exhibits post-peak strain 
softening in the macroscopic, homogenized, stress-strain response. This homogenized 
response is assumed to govern over a suitable volume of the material, appropriate 
to the microstructure of the material. Loss of positive-definiteness of the tangent 
stiffness tensor leads to a material instability which manifests as a localization of 
damage into the smallest length scale in the continuum problem [Bazant and Cedolin 
(1991)]. In GMC or HFGMC this is a single subcell. Thus, the post-peak softening 
strain energy is dissipated over the volume of the subcell that the damage localizes 
to. Since a stress-strain relationship prescribes the energy density dissipated during 
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the failure process, the total amount of energy dissipated in the subcell decreases as 
the size of the element is reduced, and in the limit zero energy is required to fail the 
structure. 

A simple way to remedy this non-physical behavior in a numerical setting is to 
judiciously scale the post-peak softening slope of the stress-strain constitutive law. 
Then, the failure energy density dissipated becomes a function of the characteristic 
length of discretized continuum. Bazant (1982); Bazant and Oh (1983) first pro- 
posed a crack band model in which post-peak softening damage (herein referred to 
as failure) in the material was assumed to occur within a band. The post-peak slope 
of the material constitutive law was scaled by the characteristic length of the finite 
element exhibiting failure; such that, the total strain energy release rate in the ele- 
ment, upon reaching a state of zero stress, and the material fracture toughness were 
coincident. In this reference, equivalence between this smeared crack approach and a 
line crack approach is presented. Subsequently, Bazant and Cedolin (1983) exhibited 
propagation of a crack band not aligned with the mesh bias. In this work, the crack 
band model is implemented within the HFGMC framework, in the MAC / GMC suite 
of micromechanics codes developed at the NASA Glenn Research Center [ Bednarcyk 
and Arnold (2002a, b)], and used to analyze crack band growth in composite RUCs. 
The following subsections provide theoretical details on the crack band model. 

3.1 Physical Behavior of Crack band 

The smeared crack band model is meant to capture the behavior of a region of a 
material wherein numerous microcracks have initiated and they coalesce to form a 
larger crack. Figure 6 displays a crack band of width w c embedded in a continuum. 
The domain of the crack band is denoted as fl' and the remaining continuum as 0. 
The crack band is oriented within the continuum such that, for a given point within 
the crack band, the unit vector normal to the crack band is n. 
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The total energy dissipated during the failure process is dissipated over fl', and the 
size w c of O' is a material property directly related to the material fracture toughness 
[Bazant and Oh (1983)]. 

_ 2Q C ( i n 

Wc crl U E t ) 

where gq is the critical stress for initiation of the post-peak regime in the ID material 
stress-strain law (given in Figure 7a), and Et is the negative tangent slope in that 
regime. The fracture toughness Qc , or critical strain energy release rate, of the 
material is given by the area under the ID traction-separation law (Figure 7b) which 
governs the cohesive response of the separation of crack faces as a crack propagates 
in the material. The energy density dissipated during failure Wp is related to the the 
material fracture toughness by the characteristic length in the material. 

Qc = w c Wf (67) 

3.2 Numerical Implementation of Crack Band Model in 
GMC/HFGMC 

The crack band model is implemented in the HFGMC micromechanical framework. 
The local subcell fields are used to govern crack brand evolution in the constituents of 
the composite. Figure 8 shows the discretization of the continuum displayed in Figure 
6. A magnified view of the crack band embedded in a single subcell is also displayed 
in Figure 8. Since the all of the energy dissipated in the crack band is smeared over 
the subcell volume, the subcell must be large enough to contain the crack band of 
width wc- Note that Figure 8 shows a 2D geometry for illustrative purposes, but the 
crack bands can also evolve in a general 3D space. 
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3.2.1 Mode I Crack Band Growth Under Tensile Principal Stress 

The orientation of the crack band in subcell by is given by the vector (see Fig- 
ure 8) and is determined from the local principal stress state d^ 77 *, v-f 1 * ) • In a 

monolithic material, cracks orient such that the crack tips are always subjected to pure 
mode I (opening mode) conditions unless there are constraints that limit the crack 
orientation. In a micromechanical analysis the composite material is composed of sep- 
arate, monolithic constituents; thus, in the matrix, the crack band runs perpendicular 
to d^ 7 \ the principal stress with the largest magnitude, > Id^ 7 ^ > jd.^ 7 ^, if 

dj^ 7 -* > 0 (tensile). Under these conditions, a crack oriented as such, is subjected to 
pure mode I loading, locally. Although, the resulting global behavior may appear to 
be mixed mode because of the influence of the fibers on the matrix crack band path. 
Crack orientation and evolution is determined differently if d^ 7 " 1 < 0 as described in 
3.2.2. The characteristic length of the subcell is determined as the dimension 
of the subcell running parallel to Crack band initiation is determined using a 

very simple, but physical, maximum stress criterion. 


-jpfi = A ^f 7) > 0 (68) 

a c 

where is the cohesive strength of the crack band. Once the crack band has 

initiated, the crack band orientation is fixed as time evolves. 

Once the crack band orientation has been calculated, the subcell compliance is 
rotated into the principal frame using the transformation matrix. 


T = [ni /37) nj 37) nj 37) ][eie 2 e 3 ] (69) 

where rif^\ and nj/* 7 '* are the principal stress directions, and ei, e 2 , and es 

are the unit basis vectors. All material degradation due to crack band evolution is 
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imposed on the rotated compliance S^\ the components of which are given by: 


0 (^ 7 ) _ rji rri rri 

u ijkl ~ ± pi- L qj^pqrs- L kr- L ls 


(70) 


The strain energy released during the formation of new surfaces corresponding 
to the growth of cracks within the crack band is assumed to be dissipated over the 
entire subcell volume. Therefore, the post-peak softening slope Ejp\ and the strain 
at which a the principal stress state is zero, is calculated using the characteristic 
length of the subcell l and the material fracture toughness (see Figure 7a). 


( 07 ) ^y_rc 

F (07)/ (07) 

a c l c 


(71) 


p(07) __ 

E IT ~ 


,(07) 


Ei\o a 


-1 


(07) 
c 


(72) 


where Eu 0 is the undamaged, axial Young’s modulus in the principal frame. It should 
be noted that e \ ^7 must be less than zero; therefore, by Equations (71) and (72), a 
restriction is placed on the maximum allowable subcell size. 


m < 26lc 7, -Biio 
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(73) 


The local, rotated, subcell strain state e\ 
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(74) 


is used to degrade the rotated compliance components. The scalar damage factor 


N AS A/TM— 20 12-217603 


20 



£>(/?7) j s calculated using the rotated strain corresponding to . 
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E ll0 A^ 


( 75 ) 


where e ^ is the value of e f' <] when the initiation criterion, Equation (68), is satisfied. 
If is less than zero, no damage occurs, and a maximum damage level of one 

corresponds to a zero stress state on the softening stress-strain curve. Also, damage 
healing is inadmissable. 

D (/?7) > 0 ( 76 ) 


Components of the rotated compliance matrix are degraded with the damage factor. 
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( 77 ) 


Since the crack band orientation is fixed upon initiation, the S^l and S^il shear 
compliances in the rotated frame are degraded, as well as the S\Ai compliance, so 
that the crack band faces normal to are free of normal and shear tractions when 
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all of the crack band energy has been dissipated (i.e. = G\c)- Once the 

compliance in the rotated frame is degraded, the compliance is transformed back to 
the global frame to yield the new subcell compliance. 


o(P 7 ) _ rji—lrri—1 o(/3^f)rri— lrp— 1 

°ijkl 1 pi 1 qj °pqrs ± kr 1 Is 


(78) 


Note that, damage introduced in the principal frame, through Equation (77), can 
induce normal-shear coupling in the global frame. 


3.2.2 Mode II Crack Band Growth Under Compressive Principal Stress 

Crack band growth under pure mode I conditions are energetically favorable; however, 
conditions arise, under which, mode I crack growth is not possible. For instance, 
the crack cannot grow under in mode I if the normal traction at the crack tip is 
compressive. This occurs when the principal stress with the maximum magnitude is 
compressive (<T| / ^ ) < 0). However, experimental data shows that monolithic materials 
subjected to global compression will eventually fracture. 

In brittle and quasi-brittle materials, it has been shown that local, internal friction 
results in shear (mode II) fracture when the monolithic material is subjected to com- 
pressive stresses [Hoek and Bieniawski (1965); Horii and Nemat-Nasser (1986); Ashby 
and Sammis (1990); Chen and Ravichandran (2000)]. Consequently, if < 0 it 
is assumed that the crack band is aligned with plane of maximum shear stress f 1 ^ 
in the matrix subcell and a Mohr-Coulomb failure criterion will be used to indicate 
crack band initiation under maximum principal compression. 


Adi) 
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< 0 


(79) 


where r^f 7 ^ is the cohesive shear strength of the matrix, and is an effective shear 


Adi) 
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stress that includes the influence of the traction normal to the crack band cr n ^\ 

f^ = \f^\4 m(Tn m (80) 

where /q is the internal friction coefficient and must lie between 0 < /q < 1.5 to ob- 
tain physically reasonable surface friction coefficient [Chen and Ravichandran (2000)]. 
Since the principal stress with the largest magnitude is compressive, the traction nor- 
mal to the maximum shear stress plane cr n ^ 7 ' ) must also be compressive. Hence, in 
Equation (80), an increase in the normal compressive stress will yield a reduction in 
the effective shear stress and an increase in apparent shear strength. Taliercio and 
Sagramoso (1995) derived relationships between internal friction coefficient, shear 
strength, the tensile strength and compressive strength of the material, yielding an 
expression for //., in terms of cr^ 7 ' 1 and Tq^. 

Similar criteria are used for shear failure under compressive loading in homoge- 
nized composite materials [ Puck and Schurmann (1998, 2002); Pinho et al. (2005)]. 
In these theories, the orientation of the crack is not aligned with the plane of maxi- 
mum shear stress because of the influence of the fibers on the crack path, but rather, 
orientation is an input to the theories obtained from experimental data. Since the 
influence of the fibers are explicitly accounted for in the present methodology, the 
crack may remain oriented with the maximum shear stress plane and the local stress 
in the matrix will drive the crack path. 

Subsequent to mode II crack band initiation via Equation (79), the compliance 
tensor of the subcell is rotated into the maximum shear stress frame using Equations 
(69) and (70), where rtf 1 '’ is a unit- vector perpendicular to the plane of maximum 
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shear stress. Features of the post-peak softening shear stress-shear strain curve, such 
as the endpoint and post-peak tangent stiffness, can be calculated from the mode II 


fracture parameters. 
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(83) 


where G '120 is the undamaged, axial shear modulus in the maximum shear stress 
frame. Again the tangent stiffness must be less than zero; therefore, Equations 

(82) and (83) place a restriction is on the maximum allowable subcell size. 
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The shear strain corresponding to the maximum shear stress ■yGi) j s obtained by 
rotating the strain tensor. 

7 (/3t) = T 1 ( f 7) e^ 7) T 2 ( f 7) (85) 

An effective shear strain 7^f 7 ^ can be defined that is work conjugate with the effective 
shear stress tEbg. 

= ( 86 ) 


where r^ 7 - 1 and d r n l ' l3y> are the shear and normal tractions acting on the crack faces 
oriented parallel to the plane of maximum shear stress, and and are the 

apparent shear and normal strains with respect to the maximum shear orientation 
in the subcell, including the effects of crack tip opening displacement. Taking the 
derivative of Equation (86) with respect to fGG yields 


df^ ] 

dfGi) 



d7 (/?7) 


(87) 


N AS A/TM— 20 12-217603 


24 



and along with the derivative of Equation (80) with respect to t 1 ' 3 " 1 ' 1 


<9rf 7) 


( 88 ) 


finally, after integrating, gives: 


(89) 


which states that the shear strain in the rotated, maximum shear stress frame and 
the effective shear strain that is work conjugate to the effective, Mohr-Coulomb shear 
stress are equivalent. 

Thus, q 1 "^ 7 ) j s US ed to degrade the rotated compliance components. The damage 
factor is calculated using 7. 


D (/?7) = 1 + 


77 - fw) 

G i2o7^ 7 ^ 


(90) 


where q^ 7 -* is the value of q^ 7 ) when the initiation criterion, Equation (79), is satis- 
fied. In the case of mode II fracture, only the shear moduli are degraded. The normal 
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direction is subjected to compression, and thus, retains its stiffness. 
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Then, the compliance tensor is rotated back to the global frame using Equation (78). 


4 Subcell Mesh Objectivity 

The key advantage of using the smeared crack band approach is that it has been veri- 
fied to provide mesh objective results within an FEM setting; this can also be said for 
implementation of this model within HFGMC. To illustrate this, a uniform, uniaxial, 
tensile displacement u was applied to monolithic, doubly-periodic, unit square RUCs 
in the .iq-dircction. As shown in Figure 9, the levels of subcell refinement ranged 
from: 35 subcells x 35 subcells, 45 subcells x 45 subcells, 55 subcells x 55 subcells, to 
65 subcells x 65 subcells. The elastic properties were chosen to match data for in- situ 
MY750/HY917/DY063 epoxy matrix [Hinton et al. (2004)]. The fracture properties 
were chosen such that the problem would demonstrate significant post-peak energy 
dissipation. All properties are given in Table 1. A single subcell, colored red in Figure 
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9, was given a 10% lower initiation strain then the rest of the domain. Additionally, 
only the subcells that lie along a horizontal line running through the weaker subcell 
were allowed to fail because as the damaged region grows, the opposite ends of the 
damaged region would interact with each other (due to periodicity), as a result of 
the periodic boundary conditions, altering the damage path and resulting in damage 
growth that is not self-similar. Restricting the damage path ensures that only subcell 
dimensions would influence the results from case to case. 

The resulting load P in the ^ 3 -direction is plotted in Figure 10 against the applied 
displacement for the four different levels of subcell refinement. It is clear that the total 
strain energy release rate dissipated in the system, and the ultimate load, is insensitive 
to the dimension of the subcell size. Of course, in more complicated RUCs, a certain 
level of subcell refinement is necessary to capture the local fields accurately; however, 
that behavior is not pathological subcell mesh dependence because the solution is 
bounded. 


5 Details of Numerical Models 

The smeared crack band model, presented in Section 3.2, is implemented within 
HFGMC, and it is used to investigated the progression of post-peak softening dam- 
age within a unidirectional fiber-reinforced composite RUC. The main objective was 
to capture the failure evolution of a composite RUC under transverse tension and com- 
pression. In tension, the lamina is susceptible to transverse cracking, as evidenced by 
the scanning electron microscope (SEM) image of a unidirectional composite failed 
under tension in Figure 11. The behavior of this mechanism is observed to be quite 
brittle, as is evidenced by tension tests on 90° coupon laminates. However, in com- 
pression a network of shear bands develop, accompanied by some matrix cracking 
and fiber-matrix debonding (see the SEM image in Figure 12) and the stress-strain 
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response if much more ductile. To emulate these failure modes, a doubly-periodic 
HFGMC RUC of the 2-3 plane of the composite was created (see Figure 13). This 
RUC was subjected to global tensile and compressive, transverse, uniaxial strains: 
±e 2 2- The constitutive behavior of the matrix subcells follows linear elasticity, cou- 
pled with the crack band model for post-peak strain softening formulated earlier. 
This model is verified against experimental data in Section 6. 

A representation of the 2-3 plane of a typical composite lamina (where the 1-axis 
runs in the longitudinal direction of the fibers and the 2- and 3-axes span the plane of 
transverse isotropy) was produced to examine the evolution of transverse cracks within 
the composite. The blueprint for the discrete model is created by randomly placing 
thirteen circular fibers with a diameter of 5 /im in a square box. The domain was 
then discretized into a number of subcells. The considered formulation of HFGMC 
only admits rectangular subcells in doubly-periodic models (parallelepiped in triply- 
periodic models); so, some of the fiber domains would overlap or lie directly adjacent 
to one another. An isoparametric formulation of HFGMC has been recently developed 
by Haj-Ali and Aboudi (2010). Mesh objective, subcell interfacial debonding has been 
formulated for GMC (see Pineda et al. (2010a, b) ) but has not been fully developed 
within the HFGMC framework. As a result, thin matrix channels were inserted 
between any adjacent fibers to avoid arrest of the crack band for non-physical reasons 
such as the inability of two adjacent fiber to separate. The final architecture and 
discretization is shown in Figure 13, which contains 81 subcells x 85 subcells. The blue 
subcells indicate fiber material and the green represent the matrix. The dimensions 
of the RUC (21.25 /im x 21.25 /irn) were chosen such that the final fiber volume 
fraction in Figure 13 would be preserved at 58%, corresponding to the experimental 
data provided in Hinton et al. (2004). 

The elastic properties for a Silenka E-glass/MY750/HY917/DY063 lamina were 
taken from data from the world wide failure exercise WWFE [ Hinton et al. (2004)]. 
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Elastic properties for Silenka E-glass fiber and MY750/HY917/DY063 epoxy matrix 
used are given in Table 2. Note that the matrix properties were correlated to represent 
the in-situ properties of the matrix, which differ significantly from those of the neat 
material [Ng et al. (2010)], and match the global properties of the composite. As a 
result of this calibration, the isotropic relationship between the shear modulus and the 
Young’s modulus and Poisson’s ratio of the matrix is not maintained. The fracture 
properties (also Table 2) were calibrated to produce and ultimate transverse tensile 
stress of 40 MPa, and the compressive properties were calibrated to match transverse 
compression data, reported in Hinton et al. (2004). 

6 Results 

6.1 Transverse Tension 

A comparison between the tensile results obtained from the HFGMC model and the 
experiment is provided in Figure 14. Only one data point is given in Hinton et al. 
(2004), implying that the tensile response is linear until ultimate failure. The first 
peak in the 022-^22 obtained from HFGMC at d 2 2 = 39.76 MPa and e 2 2 = 0.00244 
was calibrated to match the experimental data, and thus, correlates well. This is 
expected since the failure parameters were calibrated to match the experimental data. 
Subsequent to the first peak, the model exhibits a sudden drop in the global stress; 
however, the RUC then continues to reload in a nonlinear manner. A second peak is 
achieved at an applied transverse strain of 0.00308 and a transverse stress of 37.97 
MPa, which is slightly less than the stress at the first peak. After the second peak, 
the load carrying capability of the RUC is severely diminished. 

The evolution of the tensile crack band in the RUC is presented in Figure 15. At 
022 = 37.84 MPa a crack band initiates near the top of the RUC between two fibers 
(please refer back to Figure 13 for the fiber-matrix architecture). At the ultimate 
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stress 39.76 MPa, the initial crack band has grown significantly into a fully developed 
crack band that propagated across the periodic boundary into the bottom of the 
RUC. The tips of the initial crack band are arrested in lower stress, matrix rich 
regions. At <r 2 2 = 35.09 MPa, after the ultimate stress has been achieved and the 
RUC is reloading, A new crack band initiates between two adjacent fibers. When the 
second peak is reached, the subcells composing the second crack band have degraded 
significantly, and the original crack band has extended further. As the stress drops, 
Figure 15e, a third crack band initiates. The newest crack band develops rapidly, and 
the second crack band becomes stagnant. When the RUC has completely failed and 
can no longer sustain any load, the first and third crack band have bridged to form a 
nearly continuous crack that has maneuvered around the fiber inclusions and spans 
the entire height of the RUC. This failure path closely resembles the SEM image of a 
transverse crack in the glass/epoxy composite shown in Figure lib, and although the 
model does not account for fiber-matrix debonding, the model exhibits appreciable 
matrix failure in matrix subcells adjacent to fiber subcells. 

6.2 Transverse Compression 

The same HFGMC RUC is loaded in transverse compression and the stress-strain 
response is compared to the experimental data in Figure 16. The nonlinear behavior 
of the model and experiment correlate well until the model exhibits catastrophic 
failure at -104.8 MPa, which is well below the reported ultimate compressive strength 
of 145 MPa. 

Investigations into the failure progression within the RUC reveal the cause of 
the premature ultimate failure. Figure 17 shows the progression of mode II failure 
{Dc), resulting from maximum, compressive principal stresses (left column), mode I 
failure (Dt), resulting from maximum, tensile principal stress (middle column), and 
the superposition of mode I and mode II failure ( D T + Dq , right column) as the 
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RUC is loaded in global, transverse compression. The first failure initiation occurs 
between two fibers near the top of the RUC under mode II conditions at a global 
stress of -65.68 MPa. As the applied strain increases, many mode II crack bands 
form in angular matrix regions between closely packed fibers. This is exhibited in 
Figures 17b and 17c. When the global stress reaches -92.88 MPa, a mode I crack 
band initiates in the top-right corner of the RUC. At the ultimate stress, Figure 
17f both the mode I and mode II crack band have evolved. In Figure 17g, shortly 
after the ultimate stress was reached, a network of angular mode II crack bands has 
formed. These mode II crack bands are bridged by a horizontal, mode I crack band 
forming one continuous crack band that has progressed throughout the RUC. When 
the global stress has dropped severely, in Figure 17h, there is substantial matrix 
degradation resulting from mode II crack band propagation. Furthermore, one fiber 
is completely surrounded by crack bands. Comparing Figure 17h to an SEM image 
of a compressively failed carbon fiber /epoxy laminate in Figure 12b, displays closely 
resembling failure patterns, indicating that the qualitative failure mode was captured 
accurately, despite the quantitative discrepancy. 

The development of mode I crack bands corresponded to a drop in global stress 
shortly thereafter. Thus, to determine if the mode I crack bands were responsible 
for the premature failure of the RUC, the simulation was re-executed but mode I 
crack band growth was prohibited. The global < 722-022 response is presented in Figure 
16, along with the previous results and experimental data. It can be observed that 
ultimate failure is delayed considerably until the d 2 2 = -146.4 MPa, which corresponds 
well to the strength reported in Ref. [ Hinton et al. (2004)] although the ultimate 
failure strain 16% higher. Moreover, the majority of the stress-strain response of the 
RUC matches the experimental data. This indicates that, in the model, the mode II 
crack bands are responsible for the observed nonlinearity; whereas, the mode I crack 
bands are responsible for the early onset of ultimate failure, as the mode I fracture 
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toughness is an order of magnitude lower than mode II. 

Extending this observation to the behavior of actual composite would suggest 
that its non-linearity is controlled by shear band evolution, and the ultimate failure 
is due to progression of tensile cracks and fiber-matrix debonding. It may be that 
the size of the RUC (i.e. number of fibers included in the RUC) was too small to 
get an accurate representation of the quantitative response of the composite under 
compression, and the growth of the tensile crack band made an unrealistically large 
impact on the response of the RUC. Further studies evaluating the influence of fiber- 
matrix architecture on the compressive response of composites needs to be performed. 
Additionally, adjacent constraining plies with different local architectures, present in 
the experiment, may have helped to diminish the effect of formation of tensile cracks, 
which is an effect not captured in this analysis. 

The mode II crack band progression, produced by the simulation in which mode 
I crack band growth was disabled, is displayed in Figure 18. Prior to the initiation 
of the mode I crack band in Figure 17, Figures 17 and 18 are identical. However as 
the stress continues beyond the ultimate stress observed in the original simulation, 
more angular, distributed, mode II crack bands arise, as shown in Figure 18i at the 
ultimate stress of the simulation with mode I crack band preclusion. Upon ultimate 
failure, the expansive array of mode II crack bands is almost entirely adjoined and 
nearly all of the matrix subcells in the RUC have failed. The SEM image, Figure 12b, 
does not display such extensive failure. This further supports that the failure mode 
predicted with both mode I and mode II crack bands is more physically correct, and 
other factors are influencing the discrepancy between the response of the model and 
the experiment. 
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7 Conclusions 


A mesh objective, smeared crack band model was implemented at the matrix subcell 
level within the HFGMC micromechanics framework. Mode I crack bands were al- 
lowed to propagate normal to the maximum principal stress when the principal stress 
component with the maximum magnitude was tensile, and the associated strain ex- 
ceeded a critical value. If the principal stress component with the maximum magni- 
tude was compressive, however, it was assumed that mode I cracks within the crack 
band could not evolve, and instead, they grew in mode II as a result of Mohr-Coulomb 
friction, upon satisfaction of a Mohr-Coulomb failure criterion. The mode II crack 
bands were aligned with the plane of maximum shear stress. Although, the damage 
model and framework are verified using a polymer matrix composite as an example, 
the model is applicable to any brittle or quasi-brittle material including ceramics or 
concrete. 

An RUC containing 13, randomly distributed glass fibers in epoxy matrix was 
simulated under global transverse tension and compression. The results for both 
cases were compared to experimental data. The tensile response correlated extremely 
well with the test results. Although, the failure parameters were calibrated to match 
the experimental data, the successful correlation verifies that the smeared crack band 
model, utilized within HFGMC, can capture the evolution of brittle transverse cracks. 
Furthermore, the predicted failure mode mirrored SEM images of composites failed in 
tension. Under compression, the micromechanics model predicted failure well below 
the reported compressive strength. This was attributed to the formation of tensile, 
mode I crack bands in the matrix which exhibit very low fracture toughness. A 
simulation wherein mode I crack band evolution was restricted, provided quantitative 
results that more closely agreed with the experiment. However, the failure mode 
exhibited by the simulation allowing for both mode I and mode II crack bands more 
closely resembled an SEM image of a compressively failed composite compared to 
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the latter simulation. This indicates that a sensitivity study on the size and level of 
refinement of the RUC should be enacted. Future studies will determine the number 
of fibers in the RUC required to achieve convergence of the stiffness and the failure 
progression. 

Additionally, the response of the RUC to applied intralaminar shear strain 712 was 
not enacted. When subjected to in-plane shear loading numerous microcracks form 
between fibers along the fiber direction in the composite laminae [Ng et al. (2010)]. 
This 3-D effect results in the composite exhibiting a very ductile response in shear. 
To capture these 3-D, geometric effects with the proposed model, a very complex, 
3-D RUC is required. 

The RUC examples that were provided are far too computationally expensive for 
a multiscale analysis. However, insight gained from the micromechanics analysis can 
be utilized to formulate traction-separation laws, or to postulate fracture toughnesses 
that include the geometric effects and can be implemented into simpler, more com- 
putationally feasible RUCs. To preserve the stress-strain response a consistent strain 
energy density must exist across the scales. However, the mesh objectivity gained by 
utilizing the smeared crack band model at the microscale is lost if the SERR is not 
also preserved across the scales. To accomplish this, a consistent length must be used 
at both scales. Thus, the volume of the microscale RUC and the volume associated 
with the corresponding macroscale integration point must be equivalent. 
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Property 

Value 

E (GPa) 

3.7 

V 

0.35 

tc 

0.0135 

ec (weak) 

0.01215 

Qc (N/mm) 

750 


Table 1: Elastic and fracture properties used in mesh objectivity study. 
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Property 

Value 

E, (GPa) 

74.000 

v s 

0.200 

E m (GPa) 

4.65 


0.350 

Qm 

1.4 

~J W) 
e c 

0.0167 

YP 

0.0256 

Q fj ^ (N/mm) 

0.00076 

GjP'* (N/mm) 

0.00435 


Table 2: Elastic and fracture properties used in RUC simulation of 2-3 plane of a 

Silenka E-Glass/MY750/HY917/DY063 composite lamina. 
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Figure 1: Representation of triply-periodic microstructure of a composite material 

Bednarcyk et al. (2010). 



Figure 2: Discretization of a 3D, triply-periodic RUC Bednarcyk et al. (2010). 
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Figure 3: Discretization of a 2D, doubly-periodic RUC. 
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Figure 4: 2D GMC/HFGMC subcell with local coordinate frame. 



A 


Figure 5: Effects of mesh refinement on overall load-displacement response for a 

material exhibiting post-peak softening. Dashed lines indicate non-physical snapback. 
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Q 


x 5 



Figure 6: Crack band domain Q' of width w c oriented normal to vector n within a 

continuum fh 




Figure 7: Stress-strain (a) and traction-separation (b) laws governing material be- 
havior. Failure energy density Wp is related to fracture toughness Gc through the 
characteristic length. 
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Figure 8: Crack band embedded in discretized continuum. Magnified subcell displays 
crack band orientation within subcell as well as characteristic length of subcell. 



(a) 35 subcells x 35 sub-(b) 45 subcells x 45 sub- 
cells. cells. 




(c) 55 subcells x 55 sub-(d) 65 subcells x 65 sub- 
cells. cells. 




Figure 9: Four subcell meshes used in mesh objectivity study. Red subcell was given 
a lower crack band initiation strain than others. 
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u (mm) 

Figure 10: Load versus displacement for four different subcell meshes. Total strain 
energy release rate and ultimate load are unaffected by mesh refinement. 



(a) Multiple transverse cracks. [Roberts 
( 2000 )]. 


F igure 1 1 : Transverse 



(b) Magnification of a single transverse 
crack. [ Gamstedt and Sjogren (1999)]. 


cracks in glass/epoxy laminates. 
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(a) Macroscopic failure [Aragones (2007)]. (b) Localized shear bands and cracks 

[Gonzalez and Llorca (2007)]. 


Figure 12: Transverse compressive failure of carbon fiber/epoxy laminates. Macro- 
scopic failure path observed at an angle. Magnification illustrates shear banding, 
cracking and fiber-matrix debonding. 



0.005 


0.01 

X 2 (mm) 


0.015 


Figure 13: HFGMC mesh used in simulations of RUC in 2-3 plane of lamina con- 

taining 13 fibers and a fiber volume fraction of 58%. Subcells/elements occupied by 
fibers are colored blue, and subcells/elements occupied by matrix are colored green. 
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Figure 14: Transverse tensile stress versus transverse strain of E-glass/Epoxy 
posite from experiment compared to HFGMC micromechanics model. 


com- 
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(a) €22 = 0.00232, 022 = 37.84 MPa. (b) 622 = 0.00248, 022 = 39.76 MPa (First peak 

stress) . 



(c) €22 = 0.00 2 68 , 0-22 = 35.09 MPa. (d) €22 = 0.00300, 022 = 38.35 MPa (Second peak 

stress). 



(e) e 22 = 0.00 3 08, cr 22 = 37.88 MPa. (f) e 22 = 0.00316, cr 22 = 32.68 MPa. 


Figure 15: Evolution of mode I crack band in RUC subjected to applied transverse 
tensile strain. 
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x 2 (mm) 


(g) e 22 = 0.00336, (J 22 = 4.11 MPa. 


Figure 15: Evolution of mode I crack band in RUC subjected to applied transverse 
tensile strain. 



Figure 16: Compressive Transverse stress versus transverse strain of E-glass/Epoxy 
composite from experiment compared to HFGMC micromechanics model. 
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(a) 622 = -0.00403, <722 = -65.68 MPa. 



(b) 622 = -0.00490, £722 = -79.66 MPa. 



(c) e 22 = -0.00543, 0-22 = -87.71 MPa. 





(d) e 22 = -0.00578, (7 22 = -92.88 MPa. 


Figure 17: Evolution of mode I (left), mode II (middle), and superposition of mode I 
and mode II (right) crack band in RUC subjected to applied transverse compressive 
strain. 
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(e) e 22 = -0.00630, ct 22 = -99.56 MPa. 



(f) e 22 = -0.00683, a 22 = -104.78 MPa (ultimate stress). 



( g ) e22 = -0.00770, ct 22 = -96.37 MPa. 



(h) e 22 = -0.01750, ct 22 = -16.06 MPa. 


Figure 17: Evolution of mode I (left), mode II (middle), and superposition of mode I 
and mode II (right) crack band in RUC subjected to applied transverse compressive 
strain. 
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(a) e 22 = -0.00403, ct 22 = -65.68 MPa. (b) e 22 = -0.00420, ct 22 = -68.53 MPa. 



(c) e 22 = -0.00490, cr 22 = -79.66 MPa. (d) e 22 = -0.00543, cr 22 = -87.71 MPa. 



(e) e 22 = -0.00578, cr 22 = -92.88 MPa. (f) e 22 = -0.00630, cr 22 = -100.15 MPa. 


Figure 18: Evolution of mode II crack band in RUC subjected to applied transverse 
compressive strain if mode I crack band evolution is supressed. 
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(g) e22 = -0.00683, cr 22 = -106.75 MPa. (h) e 22 = -0.00770, ct 22 = -116.36 MPa. 



(i) e 22 = -0.01420, <r 22 = -146.43 MPa (ultimate (j) e 22 = -0.01750, <r 22 = -82.87 MPa. 
stress) . 


Figure 18: Evolution of mode II crack band in RUC subjected to applied transverse 
compressive strain if mode I crack band evolution is supressed. 
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